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Abstract: The three main methods currently in use for estimating the excitation class 
of planetary nebulae (PNe) central stars are compared and evaluated using 586 newly 
discovered and previously known PNe in the Large Magellanic Cloud (LMC). In order to 
achieve this we ran a series of evaluation tests using line ratios derived from de-reddened, 
flux calibrated spectra. Pronounced differences between the methods are exposed after 
comparing the distribution of objects by their derived excitation. Line ratio compar- 
isons show that each method's input parameters have a strong effect on the estimated 
excitation of a central star. Diagrams were created by comparing excitation classes with 
H/3 line fluxes. The best methods are then compared to published temperatures using 
the Zanstra method and assessed for their ability to reflect central star effective temper- 
atures and evolution. As a result we call for a clarification of the term 'excitation class' 
according to the different input parameters used. The first method, which we refer to as 
Ex„e6 relies purely on the ratios of certain key emission lines. The second method, which 
we refer to as Ex* includes modeling to create a continuous variable and, for optically 
thick PNe in the Magellanic Clouds, is designed to relate more closely to intrinsic stellar 
parameters. The third method, we refer to as E'x.[oiii]/Hf! since the [Olll]/H/3 ratio is 
used in isolation to other temperature diagnostics. Each of these methods is shown to 
have serious drawbacks when used as an indicator for central star temperature. Finally, 
we suggest a new method (Exp) for estimating excitation class incorporating both the 
[Olll]/H/3 and the He ll4686A/H/3 ratios. Although any attempt to provide accurate cen- 
tral star temperatures using the excitation class derived from nebula lines will always be 
limited, we show that this new method provides a substantial improvement over previous 
methods with better agreement to temperatures derived through the Zanstra method. 
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1 Introduction 

Much progress has been achieved in modeling the cen- 
tral stars of planetary nebulae, seen through the pro- 
gressive refinement by many authors (eg. Harm and 
Schwartzschild 1975; Schonberner 1979, 1981, 1983; 
Iben 1984; Wood and Faulkner 1986). Together with 
dynamical models for the evolution of the ejected neb- 
ula (Mathews 1966; Sofia and Hunter 1968; Kwok, 
Purton and Fitzgerald 1978; Kwok 1982), most au- 
thors are now modeling FN evolution with a more 
holistic approach (Okorokv et al. 1985; Schmidt- Voigt 
and Koppen 1987a; Stasinska 1989; 2008, Dopita and 
Meatheringham 1990, 1992; Guenther et al. 2003; 
Schwarz & Monteiro 2006; Holovatyy et al. 2008). 
This increasingly places more dependence on our un- 
derstanding of the photoionisation of the nebula and 
may affect how we interpret diagnostic tools such as ex- 
citation class which has been established to give some 
indication of the central star's effective temperature. 

Most meaningful physical parameters that can be 
derived from a PN's observed emission lines, including 
ionized and total nebular mass and the brightness and 
evolutionary state of the central star, depend on accu- 



rate distances (CiarduUo et al. 1999). This is difficult 
in our own Galaxy due to inherent problems with vari- 
able extinction and lack of central star homoegeneity 
(Terzian, 1997). The weU determined 50.6 Kpc LMC 
distance (e.g. Keller & Wood 2006), modest, 35 degree 
inclination angle and disk thickness (only ~ 500 pc, 
van der Marel & Cioni 2001), mean that LMC PNe 
are effectively co-located. Since dimming of their light 
by intervening gas and dust is low and uniform (e.g. 
Kaler & Jacoby, 1990), we can better estimate abso- 
lute nebula luminosity and size. This allows us to test 
the various methods for deriving excitation class and 
place this diagnostic on a better footing. 

We have constructed the most complete, least bi- 
ased and homogeneous census of a PNe population 
ever compiled for a single Galaxy via discoveries from 
our deep AAO/UKST Ha multi-exposure stack of the 
LMC's central 25 deg^ (Reid & Parker 2006a,b; Parker 
et al. 2005). PNe were confirmed by 2dF spectroscopy 
which gave 460 new PNe and independently recov- 
ered all 169 previously known PNe in the central area. 
These data have led to significant advances in our un- 
derstanding of the FN luminosity function (PNLF) 
(Reid & Parker, 2010) and physical parameters such 
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as temperatures, densities, nebula masses and abun- 
dances (Reid 2007). Perfect for excitation studies, the 
LMC sample provides a large range of intrinsic flux 
intensities for individual emission lines (eg. 8.8E-12 to 
3.7E-16 ergs cm^ s"^ for [O in] 5007) and interesting 
differences in line ratios. 

Despite all the advantages of using PNe in the 
Magellanic Clouds for the modeling of dynamical evo- 
lution, without deep HST imaging the central stars 
cannot be directly observed. In order to understand 
the evolution of these PNe, properties of the central 
star such as mass and luminosity must be modeled us- 
ing properties of the nebula. Fortunately, much work 
has been conducted in this area and it has been con- 
vincingly argued that the strength of emission lines in 
the nebula are uniquely determined by the properties 
of the central star (Dopita et al. 1987, 1988, 1990, 
1992). 

Much of this central star modeling depends on some 
evaluation of the so-called excitation class parameter 
which is determined from ratios of specific emission 
lines, measured from the ionised nebula. In this pa- 
per, we use the large, new population of PNe in the 
LMC to compare the three main methods currently 
used for determining excitation class. Each of these 
different methods is somewhat grounded in the physics 
of low to high ionisation states but produce differing 
results. At last, however, with a large, more repre- 
sentative population of PNe in one system, covering a 
large evolutionary range, we can observe the different 
effects produced by these differing excitation schemes. 

Having plotted, assessed and compared the exist- 
ing excitation schemes, each is found to be reflecting 
different aspects of the nebula-central star relation- 
ship, none of which is fully adequate. We suggest a new 
method for the determination of excitation class. This 
method provides a smooth transition between low and 
high excitation regimes and allows both the helium-to- 
hydrogen and the oxygen-to-hydrogen ratios to con- 
tribute a weighted proportion toward the excitation 
estimate. 

Each method is then assessed and compared by 
constructing diagrams, based on the standard H-R di- 
agram. In this case, excitation class substitutes for 
central star temperature and the H/3 luminosity of 
the nebula represents a rough estimate of stellar lu- 
minosity. The diagram constructed using the most 
correct method should provide an evolutionary map, 
over which previously determined evolutionary tracks 
should show some degree of unity if excitation class 
has any true intrinsic value. 

2 Spectroscopic Observations 
and Flux Calibration 

A major spectral confirmation program was under- 
taken in November and December 2004 comprising 5 
nights using 2dF on the Anglo- Australian Telescope 
(AAT), 7 nights using the 1.9m at South African As- 
tronomical Observatory (SAAO), 3 nights using the 
FLAMES multi-object spectrograph on the ESO Very 
Large Telescope (VLT) UT2, 7 nights using the 2.3m 



Australian National University (ANU) telescope at the 
Siding Spring Observatory (SSO) and 3 half nights us- 
ing 6dF on the UKST. In all we obtained an unprece- 
dented 7521 high and low resolution object spectra for 
LMC targets. 

Details regarding the fiuxes and fiux calibration are 
given in Reid and Parker (2010). 

2.1 Corrections for Extinction 

Extinction of light from distant objects is mainly the 

result of interstellar dust. In the case of PNe, there 
is also extinction from circumstellar dust which dims 
the light of the central star and ionised nebula. Where 
we are using and comparing line fiuxes for nebula di- 
agnostics, it is important to correct each spectrum for 
both forms of dust extinction. In the optical regime, 
the Ha/H/3 ratio was used to determine the extinction 
constant cH/3 (i.e., the logarithmic extinction at H/3 
for each nebula. For the reddening law, the interstel- 
lar extinction curves of Nandy et al. (1981) for the 
LMC were used for the extinction functions (/ (A)). 

Application of the Balmer decrement ratio 2.86 be- 
tween the Hof and ¥Ld lines was used to correct the 
spectrum in terms of the required ratios however, the 
value of c is dependent on flux calibration for each of 
these lines and any internal inconsistencies. 

The contribution of Galactic foreground dust is as- 
sumed to be low at a reddening value of E(B-V)=0.074 
(Caldwell and Coulson 1985). Full details may be 
found in Reid and Parker (2010). 

3 Determination of PN Exci- 
tation Class, E 

There are many difficulties associated with the classi- 
fication parameters for gaseous nebulae. The first de- 
termining factor is the temperature of the central star. 
The excitation class E oi & PN is a classification which 
indicates, in broad terms, the temperature of the cen- 
tral ionising star as inferred from the nebula. Ratios 
of specific emission lines provide a means of classify- 
ing PNe. High excitation PNe will be those with very 
hot stellar nuclei. Likewise, low excitation PNe will be 
those with cool stellar nuclei. Although the excitation 
class of a nebula relates closely to the central star, it 
should be viewed as a separate parameter in its own 
right which may have additional influences apart from 
the central star temperature or luminosity, including 
nebula electron temperature, electron density, mass, 
chemical composition and nebula structure. Clearly, 
for excitation, an independent, quantitative method of 
classification is required. 

The presence or absence of certain emission lines 
provide an initial indication of excitation. For exam- 
ple, the Henline at 4686Ais an initial indication of 
high excitation. In 49% (207) of the newly discovered 
PNe, this line is not present but is commonly detected 
in many of the previously known LMC PNe with only 
27% (46) non detections of He ii A4686. The presence 
of this line is an initial indication of medium to high 
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excitation of the nebula, largely ruling out possible 
confusion with Hllregions. 

Errors may be estimated and applied using a com- 
bination of the line measurement errors and the flux 
calibration errors. For our data this amounts to 10% 
for low excitation nebulae and 5% for medium- and 
high-excitation nebulae. The errors in turn affect the 
probability for any given PN to lie in a particular exci- 
tation class, but this is not true between classes where 
the presence of He ll A4686 becomes the key determi- 
nant. 

As an additional test, excitation classes were deter- 
mined for repeat observations of the same PN where 
a different 2dF field setup was used. A slight change in 
the spectrum and/or measurement of the [O in] A5007 and 
H/3 or He ll A4686 lines would be sufficient to shift any 
PN, already lying close to the edge of a class, into the 
next adjoining class. Of the 88 PNe with repeat ob- 
servations and determinations of excitation class, 63 
(71%) agreed completely with the previous class as- 
signment. Of the remaining PNe, 20 were away by 
one class, three were away by two classes (but in the 
same low-medium-high bracket) and two were sepa- 
rated further. These last 2 repeat observations suf- 
fered from lower S/N where the HellA4686 line was 
unable to be clearly measured. In each case the clear- 
est and strongest spectrum was always chosen for PN 
diagnostic analysis. 

3.1 The Exncb method 

The first classification scheme or method we examine 
for determining the excitation class of a system of PNe 
is based on the scheme originally developed by AUer 
(1956) and subsequently adopted by others (eg. Feast 
1968; Webster 1975; Morgan 1984; Gurzadyan 1988, 
1991; Acker et al. 1992). It is arguably the most 
widely used excitation scheme in the literature. Us- 
ing this method, the excitation classes are divided into 
discrete bins which may be broadly grouped into cate- 
gories described as low, middle and high. For simplic- 
ity, we refer to the two [O ill] Unes, 5007A and 4959A as 
Ni and N2 nebulae lines. 

This system is based primarily on the ratio of J(Ni 
+ N2)/-r(He II 4686A) emission line intensities for medium 
and high excitation nebula and I(Ni 4- N2)/H/3 for low 
excitation nebula where Hell is no longer detected. 
The medium-to-high ratio is sensitive to the fraction 
of O'^'^'/O'^ to be found in the nebula leading to its 
use in estimating the electron temperature Te . Where 
the He II 4686 A line is extremely weak or missing in 
low excitation PNe (class = 1-4) the line ratio I(Ni -|- 
N2)/I(Ha) is used. This ratio is linearly sensitive to 
the 0^'''/0 ratio (much like He"*" /He) and exponen- 
tially to the electron temperature Te. This is a very 
simple but effective system for excitation classification. 
We shall refer to this scheme as Ex„e6, because it relies 
solely on the measured nebula lines. 

The medium and high excitation classes are de- 
noted by integer numbers from 4 to 12-1- according to 
the log[I(Ni -I- N2)/I(Hell4686A)] ratio as shown in 
Table n 

The logarithm of the result defines the high and 



low limits of each class where the limits correspond to 
±0.1. The highest excitation PNe are placed into the 
12-|- category. They have a mean logarithmic value 
equal to 0.6. 

The medium-excitation class, _E = 4, is consid- 
ered as a transition class where the HellA4686 line 
is at the limit of detectability at 2a above the noise. 
This class may also be determined by the high ra- 
tio 7ni+N2/-^h;3 > 15 as shown in Table [T] This al- 
lows PNe to be assigned an excitation class where 
HellA4686 may be below the detectability limit but 
the 7ni+N2/-^h,9 ratio is very large. The value of the 
excitation class for any given PN is an indicator of both 
stellar temperature and PN evolution, where a radius 
estimate can be measured with a rate of expansion. 

The fractional ionization of oxygen rises and falls 
as the central star evolves. Oxygen ionization at first 
rises quickly from O^ to O^^ but then slows as it 
reaches O"''^'''. The stellar effective temperature (Te//) 
also increases with time as the average kinetic energy 
of photoionized electrons rises from 25,000 K to over 
10"" K as the star's surface temperature also rises and it 
evolves to the white dwarf phase (Dopita & Meather- 
ingham 1990; Stanghellini et al. 2002a). 

The excitation class distribution based on this Ex„e(j 
scheme is presented as a histogram in Figure [1] 

Based on the large numbers of new LMC PNe un- 
covered there are now two clear excitation classes in 
which most PNe may be found during their post AGB 
evolution. One is at excitation classes 1-2 in the low 
excitation regime and the other at classes 10-12 in the 
high range. PNe in excitation classes 1-2 have spec- 
tral types earlier than 03 on the HR diagram, which 
largely radiate photons with hv < 54.4 eV. As the exci- 
tation increases, the number of PNe decreases towards 
a trough at class 5 in the medium regime before in- 
creasing again to the high excitation peak at class 11. 
The low number density of PNe in the medium exci- 
tation regime is the major feature of this method. It 
either points to minimal change in excitation over the 
life of a PN or a very fast transition between low and 
high regimes. 

Interestingly, in the high and low excitation regime, 
the previously known PNe, which are generally brighter, 
occupy all excitation classes in somewhat similar num- 
ber ratios to the newly uncovered PNe. The exception 
to this is the medium excitation regime, where the 
brighter, previously known PNe dominate. 

The strength of the combined [O lll] line intensities 
in most of these PNe provides no indication of the 
excitation class in itself. High excitation PNe with 
Hell4686A have central stars much hotter than the 
hottest OS stars of -46,000 K (de Koter et al. 1998) 
and radiate high energy photons which produce He"*""*" 
zones which are observed as Hell recombination lines. 

3.2 The Ex, method 

An alternative approach to represent the excitation 
class as a continuous variable was introduced by Do- 
pita et al. (1990; 1992) by defining a relationship be- 
tween the ratios used in the scheme of AUer (1956). 
This resulted in a new classification scheme which is 
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Table 1: The Ex„ob Excitation- Method Classification of LMC PNe 



Low 


Middle 


High 


Class, E 




Class, E 


log(/Ni+N2/-^H/3) 


Class, E 


l0g(/Ni+N2/^H/3) 


1 


0-5 


4 


2.6 


9 


1.7 


2 


5-10 


5 


2.5 


10 


1.5 


3 


10-15 


6 


2.3 


11 


1.2 


4 


>15 


7 


2.1 


12 


0.9 






8 


1.9 


12+ 


0.6 



PNe are classified into low {E =1-4), middle {E = 4-8) and high {E ■ 
each class. The basic table layout is based on AUer (1956). 



: 9-12+) classes. The range of line ratios is shown for 



represented by: 

Sex. (0 < Sex. < 5.0) = 0.45 
Sex. (5.0 < Sex. < 10) = 5.54 



S(5007) 



S(H/3) 

f (4686) 
S(H/3) 



+ 0.78 



(1) 



(2) 



This definition again relies heavily on the strength 
of the He^'*' ion in classifying medium to high exci- 
tation PNe. For convenience, we shall refer to this 
scheme as Ex,, because it attempts to create a vari- 
able to closely represent stellar temperature. There 
are other obvious differences in the input of the two 
schemes so we compare them using the large number 
of PNe now available in the LMC. Figure [2] shows the 
Ex* scheme again separated into previously known, 
true, likely and possible PNe and binned into discrete 
classes in order to compare the scheme directly to the 
Exnob scheme shown in Figure [1] It is clear that Ex* 
moves the peak bins found in Figure [T] closer to the 
centre of the graph. It also has the effect of smoothing 
or leveling out the peaks and troughs. 



3.3 The Ex 



\oni\jnfi 



method 



As a variation on Exnob, we also measured the exci- 
tation class using (-Foiiisoo? + Soiii4959)/-FH/3 and ex- 
tended the equation across the entire luminosity range. 
This is the least used method but has been adopted 
by Stanghellini et al. (2002a) where high-excitation 
He II lines were not available. We hereafter refer to this 
method as Ex[oiii]/h/3- The graph for this method is 
shown in Figure[3]where objects containing He ll A4686 
have been plotted separately. It is clear that apart 
from the very highest and lowest bins, this method 
is generally insensitive to the presence of the high- 
excitation He II A4686 line. It does however provide a 
smooth transition across the whole evolutionary range 
of PNe based on the [Oiii]/H/3 ratio alone. As an 
indicator of central star temperature, however, this 
method should not be used. Section [4] will show that it 
is undoubtedly the least useful of the three methods. 



4 Comparison of the 3 excita- 
tion methods 

Figures 1-3 clearly show that the three methods pro- 
duce very different results. The differences are primar- 
ily due to the differing line ratios used for the medium- 
to high-excitation regimes. In all three cases, low ex- 
citation PNe are classified by [Olll]/H/3 in some form 
or other. For low excitation PNe, both Exncb and Ex* 
methods generally agree despite the use of slightly dif- 
ferent line ratios. Individual PNe fall into correspond- 
ing excitation bins with a 72% agreement. The Ex* 
method provides 5 low-excitation bins compared to 4 
for Exneb. If Bin in Ex, is combined with Bin 1, a 
direct comparison of the 4 low- excitation bins yields 
an 87% agreement. Only the Ex[oiii]/H/9 method is 
different since it permits low excitation PNe to spread 
across the full bin range without restriction. Figure [3] 
shows that He ll A4686 may be found at low excita- 
tion levels and PNe with no He ll A4686 may be found 
at high excitation levels. If we take the presence of 
He II A4686 as a strong indicator of medium-to-high ex- 
citation, then this method of excitation classification 
should be avoided. 

The biggest difference between Exneb and Ex* lies 
in the medium to high-excitation regime. The Exneb 
method keeps a consistency between low and high ex- 
citation by retaining the numerator (the [O ill] A4959 
and A5007 line fluxes) across the whole scheme. At 
medium and high excitation levels the method intro- 
duces the HellA4686 line as a denominator. This has 
been done in order to map a steady increase in the 
He+ against 0++ species. The fact that both of these 
lines are likely to increase with increasing temperature 
is problematic. Added to this are the effects of metal- 
licity on the strength of the [O ill] line (Dopita et al. 
1992). 

In the Ex, method, it is the denominator which is 
kept as a constant across the whole scheme. The H/3 
line is used as a base, against which the HeiiA4686 
line is allowed to increase, representing the increase in 
Teff of the central star. The problem with this method 
is that H itself increases with temperature and can 
even be used as a broad indicator of stellar tempera- 
ture (Dopita et al. 1992). The fiux ratio between these 
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Figure 1: Graphic representation of the excitation 
class (Exnob) for previously known and newly con- 
firmed PNe across the whole LMC survey area. Clas- 
sification is dependent upon /ni+N2/^H;3 for low exci- 
tation PNe {E = i-4) and log (iNi-t-Na/J^HolMese) for 
medium {E = 4-8) and high extinction PNe (E = 9- 
f2+). Newly calculated line fluxes are used through- 
out. The low number density of PNe in the medium 
excitation regime is the major feature of this method. 
It either points to minimal change in excitation over 
the life of a PN or a very fast transition between low 
and high regimes. 
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5.54 [F(4686)/f (4861) + 0.78] 



Excitation Class Ex. 

Figure 2: Graphic representation of the excitation 
class (Ex«) based on stellar parameters for previously 
known and newly confirmed PNe across the whole 
LMC survey area. Classification is dependent upon 
0.45[J'(5007)/F(H^)] for low to medium excitation 
PNe (class = 0-5.0) and 5.54[0.78 + F(4686)/F(H/3)] 
for medium to high excitation PNe (class = 5.0-10). 
Newly calculated line fluxes are used throughout. This 
method features a large number of PNe in the initial 
bin for each equation used (class is a half bin). The 
following bins show a steady decrease. 
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Figure 3: Graphic representation of the excitation 
class (Ex*) based on stellar parameters for previ- 
ously known and newly confirmed PNe across the 
whole LMC survey area. Classification is dependent 
upon 0.45[F(5007)/F(H^)] for low to medium excita- 
tion PNe [E = 0-5.0) and 5.54[0.78 -f F(4686)/F(H/?)] 
for PNe of medium to high excitation {E — 5.0-10). 
Newly calculated line fluxes are used throughout. This 
method features a large number of PNe in the initial 
bin for each equation used (class is a half bin) . The 
following bins show a steady decrease. 



two lines can be very narrow, making this discrimina- 
tor extremely sensitive to line measurement errors. In 
addition, multiplying the result by 5.54 may have the 
potential to increase errors by the same factor. 

Clearly, the Exncb and Ex, methods are binning 
PNe of medium and high excitation very differently 
depending on the spectral lines used, calling into ques- 
tion the usefulness of either scheme. In Figure [3] we 
compare the binning of each PNe of medium and high 
excitation in order to understand why the two meth- 
ods produce such different results. The plot shows 
that E = 5 in Ex, contains a large number of PNe 
which span up to and including excitation E = 11 
in the Exnob scheme. With each increasing Ex, ex- 
citation bin, there is an increasing agreement with 
the Exnob scheme. The reason for this is that we 
do not flnd high HellA4686/H/? levels combined with 
low [O lll]/Hen A4686 levels. This factor creates the 
curved, diagonal cutoff across the centre of the plot. 
We can therefore say that the medium excitation PNe 
in Exnob correspond to medium excitation in Ex, . We 
can also say that the highest excitation PNe in Exnob 
agree with the highest excitation PNe in Ex,. In be- 
tween, however, there is some disagreement caused by 
a large number of PNe displaying high [O ill] /He ll A4686 
levels and low HellA4686/H/? levels. 

4.1 A proposed new method for de- 
termining excitation class 

When both He ll A4686/H^ and (Ni-fN2)/He ll A4686 for 
PNe of medium and high excitation are directly com- 
pared, they produce a negative slope (see Figure [S}. 
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Figure 4: Comparison of excitation class -Eex, and 
£'Ex„^b foi' PNe of medium and higli excitation. 
Whilst there is agreement between both methods at 
the highest and lowest excitation levels, the lowest ex- 
citation levels in the Ex* method include a large num- 
ber of PNe estimated to be high excitation in the Exnob 
method. The reason for this difference is the large vari- 
ation which can exist between the [O ill] AA4959,5007 
and H/? line strengths used as the basis for the Ex* 
and Exneb methods. 



The plateau on the top-left of Figure [5] corresponds to 
cases where all the helium is doubly ionized. 

The difference between these two methods reflects 
the variation between the and [O in] A4959-I-A5007 
line strengths. The [Oiii]/If/3 ratio is not at all fixed 
but also increases to some degree with excitation. This 
serves to restrain the derived excitation where [O ill] /Hell A4686 
levels are high, thereby providing an improved method 
for estimating medium and high excitation classes. The 
most effective way to incorporate this is to change the 
input parameters so that medium and high excitation 
become: 



5.54 



+ log 



while the low excitation PNe are estimated by 
0.45 ( ^] 

\ J 



(3) 



(4) 



thereby allowing the excitation class to increase in 
step with both the He ll A4686/H/3 and the [Olll]/H/3 
ratios. It also assigns an increasing weight to the 
HellA4686 line with increasing excitation. The prod- 
uct values come from Dopita & Meatheringham (1990). 
This results in a far more robust method for estimating 
excitation class and ultimately central-star tempera- 
tures. It removes the sharp cut-off between PNe of low 
and medium excitation, based solely on the presence 
of a small measure of He ll A4686, while still allowing 
its presence to boost the excitation level. The distri- 
bution of LMC PNe into excitation classes using this 
method, which we will henceforth refer to as Exp, is 
shown in Figure (6) 
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Figure 5: A comparison of Hen A4686/H/3 and 
(Ni-|-N2)/HellA4686 for PNe of medium and high 
excitation. The negative correlation shows sufficient 
scatter to be problematic. The scatter is due to the 
variations in the [Olll]/H/3 ratios from one PN to an- 
other. The plateau on the top-left corresponds to cases 
where all the helium is doubly ionized. 



5 Excitation class comparison 
using line fluxes 

Since each scheme uses different input parameters, we 
can compare the results by plotting the line fluxes 
in relation to H/3 for each excitation scheme. This 
is an important comparative diagnostic because it re- 
veals the effects of increasing excitation on the most 
relevant emission lines. Figure [7] shows the original 
Exneb scheme with [O ill] A5007 and He ll A4686 fluxes 
plotted together with a trend line. These lines are 
sensitive to excitation and are used in the classifica- 
tion scheme. The [Nelll] A3869 line has been used by 
Morgan (1984), however it was found that the scat- 
ter was rather large. We include the trend lines for 
the He I A5876 and A6678 fluxes for comparison, as a 
ratio of these with HellA4686 could potentially de- 
fine excitation. These fluxes were measured using the 
same spectra and are published in Reid (2007). The 
He lA3889 line has been excluded as it has a metastable 
lower level and can therefore be strongly affected by 
coUisional excitation and self-absorption. It is always 
omitted when estimating He abundances in nebulae for 
the same reason. 

In Figure[8]we show the [O ill] A5007 and He ll A4686 
hne fluxes as they appear using the Ex* scheme. For 
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Figure 6: Graphic representation of the excitation 
class (i?Exp ) based on steUar parameters for previously 
known and newly confirmed PNe across the whole 
LMC survey area. Classification is dependent upon 
0.45(^5007 /^H^) for low to medium excitation PNe 

(p = 0-5.0) and5.54[log(f4959+5007/-FH/3)+^46S6/-FH/3] 

for medium to high extinction PNe (p = 5.0-13^). 
There is no longer a sudden increase in the number of 
PNe coinciding with the introduction of the Helium- 
based equation from class 5. Newly calculated line 
fiuxes are used throughout. 




2 3 4 5 6 7 8 9 10 11 12 12 + 

Excitation Class (EXneb) 

Figure 7: The position of selected high and low exci- 
tation emission line fluxes with respect to H/3 for ex- 
citation classes using the Exncb classification scheme. 
Fluxes for [O ill] A5007 and He ll A4686 have been plot- 
ted together with a polynomial trend line of best fit 
for each line. Trend lines for He I fiuxes (HelA5876 
and A6678) and [Nelll]A3869 have been plotted for 
comparison but individual fiuxes are not shown due to 
confusion. In this method, there is a general increase 
in Hen fiux in relation to H/3 with increasing excita- 
tion. This increase is not refiected in the [O ill] A5007 
flux. 



comparison we also include the trend line for the un- 
seen [Neiii] A3869, HeiA5876, and A6678 fluxes. The 
trend lines compare quite well with the Exncb method, 
despite individual objects being plotted in a different 
way. There are, however, some interesting differences 
which arise from the use of different equations. The 
He II A4686 line, which forms a gradual curve with scat- 
ter in Exncb becomes a strict exponential in Ex» . This 
forces a steadily increasing number of PNe into the 
middle excitation range. This is not surprising because 
the H/3 line was used as the denominator in the equa- 
tion. It is also used as the denominator for the y axis 
in Figure [S] thereby creating the continuous exponen- 
tial curve. The same may be said for the [O ill] A5007 
line in the low excitation area of Figure [S] When the 
results are compared to the same graph presented as 
Figure 1 in Dopita & Meatheringham (1990), the same 
trends are obvious but the scatter is increased in our 
graph due to the larger sample size and luminosity 
ranges. 

In Figure |9] we plot the same line fluxes using the 
Ex[oini/H^ scheme. This method relies on a steady 
increase in the [O ill] /H/3 ratio across the entire exci- 
tation class range. With the HellA4686 fiux ignored 
in this method, it is clear that the [O ill] /H/3 tem- 
perature indicator provides no indication as to the 
He II A4686/H/3 fiux ratio. Although the minimum lev- 
els of He II A4686 fiux per excitation class somewhat 
increase with increasing excitation, the spread is gen- 
erally even across all classes. 

Finally, Figure [TD] shows the same emission lines 
plotted using our newly suggested Exp scheme. Low 
excitation PNe without measurable He ll A4686 occupy 
classes to 5 purely based on the [O ill] /H/3 ratio. All 
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Figure 8: The position of selected emission line fiuxes 
with respect to H/3 for excitation classes using the 
Ex, classification scheme. Fluxes for [O ill] A5007 and 
He II A4686 and have been plotted together with a poly- 
nomial trend line of best fit for each line. Trend lines 
for He I fiuxes (He I A5876 and A6678) and [Ne ill] A3869 
have been plotted for comparison but individual fiuxes 
are not shown due to confusion. The increase in low 
excitation PNe is based entirely on the [O ill] 5007/H/3 
ratio, resulting in the exponential increase in flux with 
increasing excitation class. Because the medium and 
high excitation classes are derived from the Hell/H/3 
ratio alone, the same exponential curve arises from the 
position this equation is applied. This has the effect 
of forcing too many PNe into the E — 5 excitation 
category. In a similar way, the [O ill] A5007 line scat- 
ters, once the low excitation equation is dropped. This 
method also fails to incorporate the [O ill] /H/3 temper- 
ature indicator. 
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10 15 20 25 

Excitation Class (Exom/Hii) 

Figure 9: The position of selected emission line 
fluxes with respect to H/3 for excitation classes us- 
ing the Ex[oiii]/H^ classification scheme. Fluxes for 
[O III] A5007 and He ii A4686 and have been plotted to- 
gether with a polynomial trend line of best fit for each 
line. Trend lines for He I fluxes (Hel A5876 and A6678) 
and [Ne ill] A3869 have been plotted for comparison but 
individual fluxes are not shown due to confusion. This 
method relies on a steady increase in the [Olll]/H/? 
ratio across the entire excitation class range. With 
the HellA4686 flux ignored in this method, it is clear 
that the [Olll]/H/3 temperature indicator provides no 
indication as to the He ll A4686/H/3 flux ratio. 



other PNe, which include HellA4686 are allowed to 
flnd their excitation class position across the entire 
excitation range, based on both their He ll A4686/H/3 
and [Olll]/H/3 ratios. This means that they are not 
restricted to begin from class 5 upwards, but may be 
placed in a lower excitation class depending on the line 
ratios. 

The spread at each excitation level is due to the 
weighted averaging between each temperature indica- 
tor. Where [O lll] is strong and He ll A4686 is weak, 
more weight is apportioned to the [O ill] lines. Con- 
versely, where HellA4686 is strong, less weight is ap- 
portioned to the [O III] lines. Unlike the Ex* method, 
on which this method is based, there is no sudden in- 
crease in the number of PNe at class 5 where the equa- 
tions meet. Both [Olll] and Hell A4686 show a steady 
increase with rising excitation. 




1 2 3 4 5 



10 11 12 13 14 15 16 17 18 

Excitation Class (EXp) 



Figure 10: The position of selected emission line 
fluxes with respect to H/3 for excitation classes 
using the Exp classification scheme. Fluxes for 
[O III] A5007 and He ll A4686 and have been plotted to- 
gether with a polynomial trend line of best fit for each 
line. Trend lines for Hel fiuxes (HeiA5876 and A6678) 
and [Ne ill] A3869 have been plotted for comparison but 
individual fluxes are not shown due to confusion. Low 
excitation PNe without measurable He ll A4686 occupy 
the classes to 5 purely based on the [O lll]/H/3 ratio. 
All other PNe, which include Hen A4686 are allowed to 
find their excitation class position across the entire ex- 
citation range, based on both their He ll A4686/H/3 and 
[Olll]/H/3 ratios. The spread at each excitation level 
is due to the weighted averaging between each tem- 
perature indicator. There is no sudden increase in PN 
at -B = 5 where the equations meet. Both [O ill] and 
HellA4686 show a steady increase with rising excita- 
tion. 



6 Excitation class comparison 
using H/3 luminosity 

If the excitation class is efficiently representing the cen- 
tral star temperature, then it can be used to produce a 
rough type of H-R diagram. For optically thick PNe, 
the absolute H/3 flux has been suggested as a reliable 
indicator of stellar luminosity (Dopita et al. 1992). Di- 
agrams using these parameters should, if the method 
is correct, provide a rough plot of PN evolution within 
a single system. Since there is a large range in CS mass 
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Figure 11: LMC PNe plotted according to their H/3 
flux and excitation class by the Ex* scheme. 



and dynamical age within our LMC sample, the results 
appear much like scatter graphs, however, there are vi- 
tal clues which may indicate which excitation methods 
are producing the most physically realistic results. 

With the maximum luminosity conversion efficiency 
into the H/3 line derived at a stellar temperature of near 
70 000 K (Dopita et al. 1992), we should expect this 
result to be reflected in our H/3 luminosities. Stellar 
temperatures (teff) of about 70 COOK should be found 
at about excitation class 6 in the Ex, scheme and class 
8 in the Exncb scheme. According to Figure 1111 this 
class does feature some of the brightest PNe in H/3. 
There are a couple of outliers which somewhat obscure 
the result. The most notable of these is the bright PN 
between classes 7 and 8. This object is SMP32, a very 
large PN with an Ha diameter of 12 arcsec including 
the outer ionised halo (Reid & Parker 2006b). 

Our comparison conflrms that most of the bright- 
est LMC PNe do occupy the E = 3-6 excitation range. 
PNe with lower central star temperatures are less effi- 
cient at converting luminosity into their nebulae, mak- 
ing them harder to find at bright luminosity levels (Do- 
pita et al. 1992). 

At log\FH/3 fainter than -13, the number of PNe ap- 
pears to increase exponentially, in agreement with ear- 
lier predictions by Henize & Westerlund (1963) and de- 
veloped by Jacoby (1989) and CiarduUo et al. (1989). 
Clearly, most PNe, no matter their birth luminosity, 
spend the largest part of their dynamical evolution 
>ldex below the brightest possible PN in the sys- 
tem. Both Exnob (shown in Figure [12)) and Ex, agree 
that there is a deficit of PNe in the medium excitation 
range (classes 4-5) , now more apparent with the more 
representative LMC PN population available (Reid & 
Parker 2006). This new finding has a huge bearing on 
the new LMC and MASH PN samples (Parker et al. 
2006; Miszalski et al. 2008). 

The Exneb scheme in particular indicates that there 
are very few faint PNe in this medium excitation re- 
gion. One possibility for this is that the medium exci- 



Log F Hp (X4861) 

Figure 12: LMC PNe plotted according to their H/3 
flux and excitation class by the Exneb scheme. Medium 
excitation classes only have strong H/3 fluxes. This ex- 
citation classification method does not allow for evo- 
lutionary motion between excitation classes due to the 
increasing deficit of faint PNe either side of _E = 5. 



tation range is a fast, transitory period in the evolution 
of certain PNe. Not every PN may have the mass to 
enter this class but those that do will predominantly 
be young, optically thick PNe. 

The distribution of PNe at the bright end of both 
the Exnob (Figure[T2I and Ex, (Figure[TT|| plots logFna > 
— 13 are quite similar, however, there is a consider- 
able difference in excitation-class distribution for PNe 
fainter than log_FH;3 = —13.5. The difference is due 
to the use of the [O ill] A5007 line as a denominator 
in the Exncb scheme. This causes the difference to be 
more pronounced in the mid-to-high excitation regime. 
Because both the input numerator and denominator 
for objects of medium to high excitation in the Exneb 
scheme are lines which increase with excitation, they 
push plotted objects towards high excitation levels and 
prevent objects being plotted where excitation and lu- 
minosity decreases. The stronger [O ill] A5007 in com- 
parison to He II A4686, the lower the derived excitation. 
For this reason alone, we would cast suspicion over 
the use of the Exneb scheme for determining excitation 
classes as an indication of central star Tcs- 

We also produce a diagram to compare the exci- 
tation class using the Ex[oiii]/h;3 classification scheme 
where this equation is extended across the entire lumi- 
nosity range. Figure [T51 shows the result of this scheme 
compared to the H/3 fiux. Because there is no change 
of input parameter at the mid-excitation range, there 
is a commonality between the highest and lowest plot- 
ted points. This scheme therefore provides an estima- 
tion of the shape of the excitation distribution using 
a single input ratio. The decreasing number of PNe 
found with increasing excitation is more in keeping 
with what we would expect to find from photoioni- 
sation modeling. There are indications, however, that 
this method, while internally consistent, is not produc- 
ing a realistic evolutionary diagram. We would expect 



10 



Publications of tlie Astronomical Society of Australia 



^ 20 - 
X 18 - 




-16 -15 -14 -13 -12 

Log F Hp (^4861) 

Figure 13: LMC PNe plotted according to their R/3 
flux and excitation class by the Ex[oiii]/h0 classifica- 
tion scheme. In this plot, the H/3 flux is strongest 
between about excitation classes 15 and 18 but this is 
relatively equal to class 8 in Ex,. 
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Figure 14: The Ex, method for deriving excitation 
class with LMC PNe plotted according to their H/3 
flux. The diagonal lines closely represent constant ex- 
pansion velocity while evolutionary tracks show the 
computed cycle for 0.7 and 0.6 stellar mass central 
stars. From the graph, we can expect young PNe to 
brighten as their nebula shells accelerate. After the 
point of maximum luminosity, expansion begins to de- 
cline in relation to mass. 



to flnd many more PNe above ~ class 17 with fluxes 
logi^H^ < — 14 since this is where we would expect to 
find objects of high excitation which are optically thin 
and therefore have much lower H/3 luminosity. Most 
of the faintest PNe should also occupy the lowest 3-4 
excitation classes. 

Finally, we produce a diagram using the new Exp 
scheme (Figure fTS)) developed here. The brightest PNe 
(brighter than log Fuf) ~ —13) trace an upward track 
from class to class 8 before decreasing in their H/3 
flux. At that point, they either move slightly down- 
wards to join the main body of PNe or move slightly 
upwards, where PNe are increasingly optically thin. 
PNe with excitation classes greater than class 12 are 
extremely optically thin, making their true excitation 
extremely hard to plot using this method. A differ- 
ent excitation scheme should likely be applied to these 
objects. Unlike the previous excitation schemes, how- 
ever, they are now appearing in the correct part of 
the diagram. In addition to other improvements, most 
of the faintest PNe are now found in the lowest four 
excitation classes. 

PNe may be born, evolve and dissolve at various 
points within these diagrams. Bearing this in mind, 
the Exp scheme still allows us to visualise evolution- 
ary tracks where PNe move upward at the right-hand 
side of the diagram, increasing in temperature, mass 
and luminosity. After reaching their brightest poten- 
tial at the extreme right of the plot, they undergo a 
turnover where mass gain gives way to gradual mass 
loss accompanied by a steady decrease in luminosity 
and excitation class. 



6.1 A comparison with evolutionary 
tracks 



In Figure [14] we show tracks for optically thick PNe 
using the Ex, method. The tracks are derived in a 
similar fashion to that described in Dopita & Meather- 
ingham (1990). The tracks move from a density of 
50 000 cm~^ to approximately 5 000 cm~^. The 0.7 Mq 
with the highest Tctt has also been extended to where 
it would be optically thin. We include a bright exten- 
sion of the tracks in class 5, which Dopita & Meather- 
ingham (1990) consider may exist due to atmospheric 
extinction of the He^ ion in this particular tempera- 
ture range. Although we have included it, we don't see 
any real evidence for it in the Ex, version of the excita- 
tion scheme. Using the Ex[oiii]/h/3 method, however, 
shown in Figure [131 SMP32, the very bright object sit- 
ting by itself in the other two plots may now be seen 
as part of an evolutionary loop or bright extension. 

In Figure [TS] we show the diagram using the Exp 
method in order to overlay the same evolutionary model 
shown in Figure [14] using the Ex, method. Estimated 
mass ranges are extended over 13 excitation classes 
rather than 10 (see Figure [T4)l due to the different form 
for deriving excitation class. Extreme optically thin 
PNe, above excitation class 12 and < magnitude -14, 
have not been plotted as their true positions cannot 
be correctly plotted using current methods. The po- 
sitions of PNe on this plot correspond very well with 
the positions predicted using the evolutionary model 
(Dopita et al. 1990). 
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Figure 15: The Exp method for deriving the excita- 
tion class with LMC PNe plotted according to their 
H/? flux. The diagonal lines closely represent constant 
expansion velocity while evolutionary tracks show the 
computed cycle for 0.7 and 0.6 stellar mass central 
stars. Estimated mass ranges are extended over 13 
rather than 10 (see Figure I14p excitation classes due 
to the different form for deriving excitation class. 



Figure 17: Temperatures derived using the Zanstra 
method (Villaver et al. 2003; 2007) are plotted against 
excitation classes for the same CSPNs derived using 
the Ex„e6 and Ex[oiii]/Hp methods. A correlation co- 
efficient of 0.63 for Ex„ei, and 0.37 for Ex[oiii]/h/3 is 
found. Data for error bars are from Villaver et al. 
(2003; 2007). 
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Figure 16: Temperatures derived using the Zanstra 
method (Villaver et al. 2003; 2007) are plotted against 
temperatures for the same CSPNs derived using the 
Exp and Ex* methods. A correlation coefficient of 
0.742 for Exp and 0.616 for Ex, is found, although 
Ex-based methods luainly return teiuperatures higher 
by a mean of 26,000 ±31,000 K for Exp. The fine repre- 
sents the 1:1 line of unity between both methods. Data 
for error bars are from Villaver et al. (2003; 2007). 
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Figure 18: Temperatures derived using the Zanstra 
method (Villaver et al. 2003; 2007) are plotted against 
excitation classes for the same CSPNs derived us- 
ing the newly established Exp method. We show a 
separation between optically thick and optically thin 
PNe where the [Oll]3727/H^ and [Nn]6584/Ha are 
below unity at orders less than 1.5 (Kaler & Ja- 
coby, 1989; Jacoby & Kaler, 1989). In addition, 
Te//Z(Hen)/Te//Z(H/3) > 2 for these PNe also in- 
dicates an optically thin nebula. A correlation coef- 
ficient of 0.924 for Exp is found using optically thick 
PNe. The line represents the 1:1 line of unity between 
both methods. Data for error bars are from Villaver 
et al. (2003; 2007). 
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7 A comparison with published 
temperatures 

PN central star (CSPN) temperatures estimated through 
the excitation class method are compared to the most 
reliable CSPN temperatures for the same objects in the 
published literature from Villaver et al. (2003; 2007). 
Altogether these papers report temperature estimates 
for 65 CSPNs in the LMC. A total of 45 of these ob- 
jects are within our survey area and are therefore used 
for comparison. 

To gain CSPN estimates, ViUaver et al. (2003,2007) 
have used the method originally developed by Zanstra 
(1931). Although this method has been widely used 
(Harman & Seaton 1966; Kaler 1983; Stanghellini et al. 
2002b), the use of different nebula recombination lines 
produces vastly different temperature estimates. For 
example, the Zanstra temperature for LMC PN J33 is 
38 700 ±5 400 K using the H/3 line and 86 700±5 400 K 
using the He ii A4686 line. The difference between the 
two temperatures is the well-known 'Zanstra discrepency' 
effect (Kaler 1983; Kaler & Jacoby 1989). Our compar- 
ison tests using Zanstra temperatures calculated with 
the H/3 line, where the nebula is known to contain 
a HellA4686 line, show scatter plots with no corre- 
lation for temperatures above 60 000 K. Temperatures 
derived using the He Il4686 line are undoubtedly more 
reliable as nebulae are optically thick to He ll ionising 
photons (Villaver et al. 2007). 

Whereas our spectroscopic data are homogenous, 
this is not quite true of the data used in Villaver et al. 
(2003,2007). Stellar photometry was obtained from 
HST using the Space Telescope Magnitude (STMAG) 
system. Extinction constants were not available in the 
literature for all the PNe; 6 H/3 fluxes were from differ- 
ent authors and all He ll A5686 fluxes were gained from 
5 different authors. Nonetheless, with careful treat- 
ment, the resulting temperatures probably represent 
the most accurate currently available in the literature 
to date. 

In Figure [T6l we show a comparison of CSPN Te//s 
using the Exp and Ex» methods against the T^z/s us- 
ing the Zanstra method as published by Villaver et al. 
(2003; 2007). Transformation from excitation class to 
temperature was made using; 



LogTe// = 4.439 + [0.1174(Exp)] - [0.00172(Exp)^ 



(5) 



which is based on the transformation given in Do- 
pita et al. (1992) but adjusted to match the published 
Zanstra temperatures. 

The published Zanstra-based temperature of 372 400±97 
(HeiiA4686) or 526 000±158 500 K (H/3) for SMP93 
have been omitted as our temperature for this object is 
estimated at 96 600 K (Exp) or 147 600 K (Ex.). Like- 
wise the published, hydrogen-based Tcff Z of 17 100±1 100 K 
for SMP28 has been omitted since we find a strong 
HellA4686 line present, raising the temperature to 
81 800 K according to the Exp result. 

Figure [TB] shows a modest correlation between tem- 
peratures derived using the Zanstra method and those 
derived using the Exp and Ex» excitation class meth- 
ods. On average. Ex temperatures derived using Exp 



are 26 000 ±31 000 K higher than those for the same 
CSPNs using the Zanstra method. Much of which is 
due to 9 CSPNs between 50 000 and 100 000 K which 
have ToffZ between 3 and 0.5 times lower than indi- 
cated by the Exp method. This phenomenon has previ- 
ously been observed in Zanstra-based temperature re- 
sults (Kaler & Jacoby 1989; Stasinska & Tylenda 1986; 
Schonberner & Tylenda 1990; Gruenwald & Viegas 
2000). It is also analogous to the test of Zanstra tem- 
peratures conducted by Mendez, Kudritzki & Herrero 
(1992) where Zanstra temperatures were compared to 
spectroscopic Tcff. Low Zanstra temperatures were 
found for 6 CSPNs in the 60 OOOK to 90 OOOK range. 
In each case the Zanstra temperature approaches unity 
as the temperature increases. 

In Figure[l7]we show a comparison of temperatures 
using the Zanstra method with excitation classes de- 
termined using both the Exnob and Ex[oiii]/h/3 meth- 
ods. With both methods, there is a gap between PNe 
of low and high excitation. Low excitation PNe occupy 
Tcff Z between 20 000 and 80 000 K while extending up 
to excitation class 9. High-excitation PNe occur from 
class 4 upwards however occupy almost any tempera- 
ture at any given excitation class. No linear fit between 
excitation class and temperature is possible, however, 
a logarithmic fit with errors up to 100% may be pos- 
sible. 

Finally, we plot our new Exp method and sepa- 
rate PNe which are deemed to be optically thin due to 
a mean [Oii]3727/H/3 = 0.2±0.08 and [Nii]6584/Hq 
= 0.3±0.23 with the [O ll]3727/H/3 ratio providing the 
stronger indicator fFigure [T8)l . Plots shown by Kaler & 
Jacoby (1989) and Jacoby & Kaler (1989) indicate any 
N/0 ratio below 0.26, corresponding to [Nll]6584/Ha 
= 0.20 will indicate an optically thin nebula. They 
also allow for a factor of 5 enrichment in nitrogen by 
selecting N/Ha = 1 as a secondary criterion (Kaler & 
Jacoby 1989). Another indicator of optically thin PNe, 
reffZ(Hell)/r<,ffZ(H/3) > 2 (e.g. Kaler & Jacoby 1989), 
is also in agreement for these PNe. A mean agreement 
of 83±22% between Zanstra and Exp temperatures and 
a correlation coefficient of 0.924 for Exp are found us- 
ing only optically thick PNe. This shows that, as long 
as optically thin PNe can be distinguished, the Exp 
method will produce very good temperature estimates. 
Figure [TH] gives us confidence that the Toff of a CSPN 
can be determined from emission line ratios in opti- 
cally thick PNe with some accuracy using the new Exp 
method. Arbitrarily, the graph indicates that temper- 
atures estimated by excitation class for optically thin 
PNe should be reduced by at least 50%. 

600 K 



8 Conclusion 

The utility of the commonly used excitation class es- 
timator for PNe is investigated and appraised using a 
large, new sample of PNe in the LMC. 

The importance of determining the best form for 
deriving the excitation class is going to be an impor- 
tant part of our on-going study of LMC PNe. The large 
number of PNe now available over a wide luminosity 
and evolutionary range has permitted the excitation 
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class parameter to be re-investigated and re-worked. 
The 3 main methods current in the Uterature for de- 
termining excitation class have been carefully plotted, 
compared and evaluated. We advise against using the 
Exnob scheme in its current form due to the distor- 
tion it creates at the medium-to-high excitation region. 
The Ex* scheme also appears to create some anomalies 
at the medium excitation levels where the input equa- 
tion is re- worked. Using ([Olll] A4959 + A5007)/H/3 for 
all the PNe produces a smooth transition but suffers 
for not applying the HellA4686 ratio. 

The new Exp method offered in this paper ap- 
pears to improve excitation and stellar temperature es- 
timates by allowing both He ll A4686 and [O ill] A5007 
to become constant variables against the H/3 intensity. 
It also minimises any sudden increase in PNe at the 
point where the Hen line is introduced into the de- 
rived equation for high-excitation PNe. The method 
also correlates well with the best published Zanstra 
temperatures, placing this new excitation class param- 
eter on a far more solid footing. 

This presentation has been prepared in order to lay 
the groundwork for further research using central star 
temperatures, electron densities, nebula mass and ex- 
pansion velocities and new photoionisation codes. The 
disparity between the current methods for assigning 
excitation classes has been highlighted in this paper. 
Our suggested Exp method appears to produce encour- 
aging results as it allows excitation to increase using 
both the [Oiii]/H^ and HeiiA4686/H^ OS tempera- 
ture indicators. 

It is clear that excitation is a quality that remains 
key to the ionised nebulae well into the ageing process 
and must be correctly modeled. Guided by these re- 
sults and comparisons, we will continue to model the 
excitation class based on the ratios of the same ionic 
species, for example the ratios of O"*" to O^"*" and He"*" 
to He"*""*". This method, used already by some authors 
(e.g. Martin & Viegas 2002) is expected to be a good 
CS temperature indicator. Our preUminary results us- 
ing the LMC PNe are encouraging. With new pho- 
toionisation models, we intend refining the input pa- 
rameters and anticipate that new evolutionary tracks 
and diagrams describing the AGB to WD stages will 
be produced. 
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